Heterogeneous Dynamics of Coarsening Systems 
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We show by means of experiments, theory and simulations, that the slow dynamics of coarsen- 
ing systems displays dynamic heterogeneity similar to that observed in glass-forming systems. We 
measure dynamic heterogeneity via novel multi-point functions which quantify the emergence of dy- 
namic, as opposed to static, correlations of fluctuations. Experiments are performed on a coarsening 
foam using Time Resolved Correlation, a recently introduced light scattering method. Theoretically 
we study the Ising model, and present exact results in one dimension, and numerical results in two 
dimensions. For all systems the same dynamic scaling of fluctuations with domain size is observed. 



Glassy and jammed materials display similar phe- 
nomenology, characterized in particular by slow and 
nonequilibrium dynamics, whose microscopic origin is 
still being actively investigated 0, 0]. Recent re- 
search has shown that despite the absence of static order- 
ing, glass-formers exhibit non-trivial spatial correlations 
of the local dynamics, resulting in dynamic heterogene- 
ity 0, ■ Here, we take the view that slow dynamics 
is intrinsically associated to dynamic heterogeneity, as 
suggested by studies of the glass transition [fj. 

Dynamic heterogeneity has to be measured by means of 
statistical correlators that probe more than two points in 
space and time. An ideal experiment or calculation would 
compare local configurations around position r at times 
t and t + At via a two-time quantity, F(r, t, At). Tradi- 
tionally, only the dynamics averaged over t, r, or thermal 
histories is discussed, F av (At) = J v d d rF{r,t, At)). 
In a glass-former, F av could be for example the self- 
intermediate scattering function at a given wavevector. 
By contrast, our goal is to detect spatial correlations of 
the local dynamics. A natural correlator is 



C(r,At) 



1 

V 



d d r'(F(r',t, At)F(r'+r, t, At)}-** (At), 



(1) 

built from two-point, two-time quantities. It is easier to 
measure the volume integral of JU, a dynamic suscepti- 
bility x(At) = J v d d rC(r, At), which can be rewritten as 
the variance of the fluctuations of the two-time dynamics, 



x (At) = V 




d d rF(r, t, At) 



Fl(At) 



(2) 

Physically, dynamic fluctuations increase when the num- 
ber of independent dynamic objects decreases, but nor- 
malizations ensure that \ remains finite in the thermo- 
dynamic limit, except at a dynamic critical point, as dis- 
cussed for supercooled liquids Q. 



In this paper, we focus on coarsening phenomena, as 
one of the simplest and better understood physical situ- 
ations characterized by slow dynamics and aging Q . We 
introduce new experimental and theoretical methods to 
access x an d C, and show that coarsening systems dis- 
play heterogeneous dynamics similar to that observed at 
equilibrium in glass formers. We perform experiments on 
a coarsening soft material, a dry foam, and we study the 
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FIG. 1: Dynamic susceptibility x(tw,At) measured experi- 
mentally in a coarsening foam (top), and numerically in the 
d — 2 Ising model (bottom) for various ages t v . 
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Ising model in one and two dimensions. Although very 
different, these systems display the same dynamic scaling 
of fluctuations with domain size. 

Experimentally, soft materials are well suited for in- 
vestigating dynamic heterogeneity because the relevant 
time and length scales are much larger than in molec- 
ular systems, which greatly simplifies detailed dynamic 
measurements. We use dynamic light scattering in the 
strongly multiple scattering limit (Diffusing Wave Spec- 
troscopy, DWS ,9j) to probe the dynamics of a shav- 
ing cream. The foam coarsens as the time f w since its 
preparation increases. The time-averaged intensity auto- 
correlation function measured by DWS, g2(t V! ,Ai) — 1, 
decays exponentially, with a characteristic rate, r(t w ), 
that decreases with t w . This results from spatially and 
temporally localized random rearrangements |l0|. 

Traditional light scattering measurements use a point 
detector and require an extended time average of the 
intensity correlation function. Consequently, no infor- 
mation on the fluctuations of the dynamics is accessi- 
ble. To overcome this limitation, we use the recently 
introduced Time Resolved Correlation technique 01 • A 
charge-coupled device (CCD) camera is used to record, at 
a constant rate, the speckle pattern of the light scattered 
by the foam [12J. The de gree of correlation, cj, between 
speckles at times t w and t w + At is calculated as 



d(t w ,At) = 



(I P (t w )I p (t w + At)) p 
I P (t w )) p (I p (t w + At)}j 



-1 , 



(3) 



where I p (t w ) is the intensity measured at time t w by the 
p-th CCD pixel and (• • -) p denotes an average over all pix- 
els. Since different configurations yield different speckle 
patterns, c/(t w , At) quantifies the degree of correlation 
of the foam between times t w and t w + At. Because the 
experiments are performed in the transmission geometry 
and in the strong multiple scattering regime, each CCD 
pixel collects light coming from the whole scattering vol- 
ume. Therefore, c/(£ w ,At) provides time-resolved but 
spatially-integrated dynamical information. The stan- 
dard intensity correlation function is obtained via a fur- 
ther average, g2(t w ,At) — 1 = (c/(t w , At))y, which is 
analogous to the correlator F av defined above. Here, 
(•••)t denotes an average over [t w , t w +T], with T ^> 
but short enough to prevent any significant change of 
the dynamics due to coarsening. In pratice, we took 
T < 0.05t w for all ages, and checked that our results 
do not depend on this choice. Dynamic fluctuations are 
quantified via the variance of c/, similarly to Eq. J2J: 
X {U,At) = (cj(U,At)) T - (ci(£ w ,Ai))|. 

The dynamic susceptibility x(^w. At) measured during 
the aging of the foam, from t w — 3350 to 26200 sec is 
shown in Fig.^ Trivial contributions to the fluctuations 
due to the CCD noise and the finite number of speckles 
have been subtracted from the data For all t w , we 
find that \ exhibits a peak at time lags At*(t w ) close 
to r _1 (t w ). Moreover, as the foam ages and coarsens, 



the height of the peak, X*(*w) = x(tw,At*) increases 
and its position shifts to larger times, in striking analq 
with numerical observations in supercooled liquids ["" 
To our knowledge, no experimental measurement of the 
dynamic susceptibility was so far reported. 

A similar behavior can be observed in a very differ- 
ent coarsening system. We study the dynamics of the 
Ising model on a regular lattice quenched from a random 
state to the ferromagnetic phase. Domains of positive 
and negative magnetization develop and grow with t w , 
in analogy with the bubbles of the foam. The Hamil- 
tonian is H = — J^V (JjCTj, the sum being over nearest 
neighbor pairs. From l[T|l. one would naively study Cq(1 — 
k,t w ,At) = {a k (t)a k (t w )ai(t)ai(t w )) - {a k (t)a k (t w )) 2 , 
where t = At + t w . However, Co is trivially dominated by 
equal-time two-point correlations: consider for instance 
the large time limit, where Cq(1 — fc,i w ,At — > oo) — > 
(a k (t ) 07 (t ) ) (a k (t w ) a i (t w ) ) ■ The appropriate correlator to 
consider is instead 11 1 



C(l-k,t w ,At) 



- (cr fe (t)a k (t w )ai (t)ai (t w )) 
+ (o-fc (t)a k (t w )) (ai (t)ai (t w )) 
+ (ak (t)m (t)) (<7fc (t w )a t (t w )) 

-fo(t)<7j(tw))<0*(*w>j(*)>. (4) 



where the relative signs are consequences of the fermionic 
nature of fluctuations in the Ising model. The corre- 
sponding susceptibility is defined as in Eq. (J2J, 



x(t w ,At) = ^C(M W ,A*). 



(5) 



Note that © is not positive definite, but to ease the com- 
parison with previous research, definitions J3J and (j^J are 
chosen so that these quantities are eventually positive. 

Dynamic fluctuations in the d = 2 Ising model are 
measured in Monte Carlo (MC) simulations. The sim- 
ulated system, L, must be large enough that the mean 
domain size, R(t w ), satisfies R(t w ) <C L at all times, 
but small enough not to average out the fluctuations. 
We used L — 600 and averaged the results over 2 • 10 3 
independent initial conditions, the total simulated time 
being 10 4 MC steps. Due to computational limitations, 
we could only measure x(t W5 At) as the linear integral of 
Q measured along the x and y axes. As seen in Fig. ^ 
the observed dynamical fluctuations are similar to those 
observed for the foam. 

This similarity lies in the common physical mechanism 
responsible for the slow, heterogeneous dynamics: the 
domain growth driven by the reduction of interfacial en- 
ergy. As the system coarsens, the size of the regions that 
undergo correlated rearrangements increases, so that the 
number of independent regions in the probed volume, 
N(t w ), decreases, which in turn increases the amplitude 
of the dynamic fluctuations. One would therefore expect 
that x*(«w) ~ N- 1 ^) ~ R d (t w ). 
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To substantiate this prediction, we first consider an 
analytically solvable case. We study the non-equilibrium 
dynamics of the d = 1 Ising model. The spi n chain 
evolves according to standard Glauber rules [l5j. From 
general expressions for two-time, multispin correlation 
functions following a quench derived in [lfj , we compute 
exactly the quantities (@J and (J5J. Although straightfor- 
ward in principle, actual calculations require technically 
involved algebraic manipulations, detailed in 

Consider first equilibrium dynamics which we ex- 
pect to be dynamically homogeneous. Indeed, we find 
that the correlation |0J vanishes exactly, C cq (n, At) = 
lim^^oo C(n,t w , At) = 0, for arbitrary distances n, 
times At, and temperatures T > 0. This implies that 
two-time multispin correlations factorize into two-spin 
static correlations, and shows that the correlator (0J is 
well suited to revealing the existence of non-trivial, dy- 
namic, out of equilibrium correlations. 

The situation is more interesting in the coarsening 
regime. We focus on the scaling behavior for large times 
At, t w — ► oo and distances n — * oo, with scaling variables 
a = At/t w and 77 = n/R(t w ) fixed, where R(t w ) ~ \/tZ- 
One finds scale invariance, 

/At n \ 

where the scaling function fc(&, rj) is given in |14|. From 
©, the susceptibility JSJ) scales as 

X (t w ,At)~R(t w )Fc(j^J, (7) 

with Fc(a) — J drj fc(ot,rf). Both scaling functions are 
displayed in Fig. [2 In particular, x(£ w , At) qualitatively 
resembles the results in Fig. Q] and Eq. JJJ shows that 
the peak height scales as expected, x*(t w ) ~ R(t w ). 

To develop a more direct understanding of 10, we 
note that the dynamics of the Ising spin chain can be 
mapped to a diffusion-limited annihilation process, where 
one studies the dynamics of the walls separating domains 
of opposite magnetization rather than that of the spins 
themselves 01 ■ Trajectories of the walls are random 
walks that annihilate whenever they meet. In a space- 
time diagram 6], the spins Uk(t w ), &i(tw)i &k(t) and 
o\ (t) occupy the corners of a rectangle of size n — k — I 
by At = t — t w , and the spin products in (@J are de- 
termined by the parity of the number of random walks 
crossing the relevant edge. Labeling edges in the order 
left - right - top - bottom we denote, for example, an 
odd number of random walkers crossing the left and bot- 
tom edges as 1001. Because walls annihilate in pairs, 
the number of walkers crossing the rectangle is even, so 
that there are only 8 possible situations. In terms of the 
corresponding probabilities, Q may be reexpressed as 
C(n,t w ,At) = 8(p ioiPioio +P0110P1001 -PooooPim - 
Pooii Pnoo)- These probabilities can be evaluated to lead- 
ing order via standard random walk arguments |l4j | , when 



10~ 3 



- 1 ii, 

pier 1 


-1 1 — r| 1 1 — T| 1 1 — T| 1 1 — T| 1 1 — n 


1 1 I I I I ■ 1 1 1 1 1 I I Ml| 1 1 1 1 1 1 II 




11? 10- 4 
- Jt 

10- 8 






'-_ 10- 10 
, ,1 


1 1 . . .liC 




10" 1 rj 10 1 

1 ill 1 ill 1 ill 1 ill 1 11 



1CT 3 1Q- 1 10 1 10 3 



a 

FIG. 2: Exact scaling function of Eq. J7J for the dynamic sus- 
ceptibility of the Ising chain. Dashed lines show the asymp- 
totic behavior, Fc(a —* 0) ~ a and Fc(a — ► 00) ~ a~ 2 . In- 
set: Exact scaling function of Eq. ©, for a = 1CP 2 , 1 and 10 
(from top to bottom on the left). Dashed lines represent the 
asymptotic behavior estimated via random walk arguments. 

the number of trajectories crossing the rectangle is small, 
i.e. for diluted walls, n -C R(t w ), and short time de- 
lays, At < t w . For At < n 2 < i w , C{n,t w ,At) con- 
verges to an n-independent plateau of height (2/-7T 2 )a, 
see Fig. |21 When n 2 « A( < i w , on the other hand, 
C{n,t w ,At) w 7r- 1 (l-2/7r)n 2 /t w , implying that C grows 
like rj 2 with a Ai-independent amplitude, see Fig. [3 The 
random walk picture becomes too complicated when ei- 
ther n 2 or At are large compared to t w and we refer to our 
exact results in this regime ^i). For n 2 <C t w -c At, the 
?y 2 -dependence found above for n 2 -C At -C t w persists, 
but now with an amplitude that decreases as oT 2 . For 
large n 2 3> t w , finally, we find a Gaussian cutoff in n with 
a width of order y/tZ for both At <C t w and At ^ t w . 

The asymptotic behavior of x(*w, At) follows from the 
above discussion. For At -C t w , the integral is dominated 
by the plateau region of fc- Since the plateau grows as 
a, so does @. For At 3> t w , on the other hand, the 
decrease in the amplitude of the 7y 2 -part of C controls the 
variation of x which decays therefore as a~ 2 . Rigorous 
analysis shows indeed Fc(a — > 0) ~ 4 (V2 — l) 7r _3 / 2 a 
and Fc(a — > 00) ~ (8/5) (8^/2 - 9) 7r- 3 / 2 a~ 2 , see Fig.0 

For the d = 2 Ising model and the foam no analyti- 
cal results are available. Instead, we test the scaling of 
x(t w , At) by measuring R(t w ) and the characteristic re- 
laxation time. In simulations, R(t w ) is determined from 
the decay at large distance of the equal-time structure 
factor, while the characteristic decay rate scales as l/t w , 
r(t w ) ~ 0- I n t ne simulations, x is a linear inte- 
gral over distance, as mentioned above. Thus, one ex- 
pects x*(^w) ~ R(t w )- For the foam, R(t w ) scales as 
the bubble size, which is also proportional to the photon 
transport mean free path, £*(t w ) 0, which we obtain 
independently from transmission measurements 0- We 
expect therefore X*(^w) ~ i* 3 (t w )- 

As seen in Fig.|3| not only does the peak of the dynamic 
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FIG. 3: Dynamic scaling of all dynamic susceptibilities in 
Fig. (same symbols). Top: experiments on a foam, the 
line is the scaling function of Eq. (JHJ. Bottom: numerical 
simulations of the d — 2 Ising model. 



susceptibility scale as expected, but data for all times col- 
lapse onto a master curve, a scaling behavior analogous 
to that found for the Ising chain, Eq. (JJJ. For the foam, 
the shape of the master curve can be explained by not- 
ing that at all times, cj(t w ,At) ~ exp[— j(t w )At], where 
7(i w ) is a decay rate fluctuating with the number of rear- 
rangements, of mean T(t w ) = (7(£ w ))t- To leading order 



in the variance of j(t w ), <r 2 (t w ) = 
£* 3 (t w )r 2 (£ w ), one can estimate 



<7(*w) 5 



-<7(*w)>! 



X (U,At)~t 3 (t w )[r(t w )At] e 



-2r(t w )At 



(8) 



The solid line in the top panel of Fig. [3] shows that the 
scaling function f(x) = x 2 exp(— 2x), suggested by JSJl, is 
in very good agreement with the experimental data. 

In summary, we have defined and analyzed dy- 
namic spatial correlators and susceptibilities revealing 
the growth with time of dynamic heterogeneity in coars- 
ening systems. The increasing dynamical fluctuations 
are consequences of a reduced number of independent 
dynamical domains — bubbles or magnetized domains — 
as coarsening proceeds. The non-fractal morphology of 
the domains implies a particularly simple scaling of spa- 
tial correlators, and therefore of their volume integral, 



X(t w ,At) ~ i? d (£ w )/(r(£ w )At). However, dynamic het- 
erogeneity is also found in systems where dynamic do- 
mains are believed to have a fractal morphology 0, 0] • 
In this case, a more complicated dynamic scaling can be 
expected 0, Q ■ The techniques presented in this work 
will help identify and characterize dynamic heterogene- 
ity, an endeavor that appears as a crucial step in gaining 
a better understanding of the dynamic slowing down of 
many glassy and jammed materials. 
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